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Abstract. - We study the directed polymer with fixed endpoints near an absorbing wall, in 
the continuum and in presence of disorder, equivalent to the KPZ equation on the half space 
with droplet initial conditions. From a Bethe Ansatz solution of the equivalent attractive boson 
model we obtain the exact expression for the free energy distribution at all times. It converges at 
large time to the Tracy Widom distribution F4 of the Gaussian Symplectic Ensemble (GSE). We 
compare our results with numerical simulations of the lattice directed polymer, both at zero and 
high temperature. 



Much progress was achieved recently in finding exact 
solutions in one dimension for noisy growth models in 
the Kardar-Parisi- Zhang (KPZ) universality class [Tj,[2] , 
and for the closely related equilibrium statistical mechan- 
ics problem of the directed polymer (DP) in presence of 
quenched disorder [jjj. The KPZ class has been explored 
in several recent experiments [4j|5], and the DP has found 
applications ranging from biophysics [6] to describing the 
glass phase of pinned vortex lines [7] and magnetic walls 
(8). The height of the growing interface, h(x,t), corre- 
sponds to the free energy of a DP of length t starting at 
point x, under a mapping which is exact in the continuum 
(Cole-Hopf), as well as for some discrete realizations. Not 
only the scaling exponents h ~ i 1 / 3 , x ~ t 2 ^ 3 are known 
9.10], but also the one-point (and in some cases the many- 
point) probability distribution (PDF) of the height have 
been obtained 
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Their dependence in the initial 
condition was found to exhibit remarkable universality at 
large time, with only a few subclasses, most being related 



to Tracy Widom (TW) distributions 15 of largest eigen- 
values of random matrices. Most of these subclasses were 
initially discovered in a discrete growth model (the PNG 
model) [TT -13 which can be mapped onto the statistics 
of random permutations fll] , and a zero temperature lat- 
tice DP model 
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Recently, exact solutions have been 
obtained directly in the continuum at arbitrary time t, for 
the droplet 18 -21 , flat 22 23 and stationary [24] initial 



conditions. The PDF of the height h(x,t) converges at 
large time to F 2 , the Gaussian unitary ensemble (GUE), 



and to F\, the Gaussian orthogonal ensemble (GOE) uni- 
versal TW distributions, for droplet and flat initial condi- 
tions respectively. One useful method which led to these 
solutions introduces n replica and maps the DP problem 
to the Lieb Liniger model, i.e. the quantum mechanics of 
n bosons with mutual delta- function attraction, a model 
which can be solved using the Bethe Ansatz. 

The KPZ equation on the half line x > 0, equivalently 
a DP in presence of a wall, is also of great interest. In the 
statistical mechanics context constrained fluctuations are 
important for the study of fluctuation-induced (Casimir) 
forces 
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and for extreme value statistics. In the sur- 
face growth context one can study an interface pinned at 
a point, or an average growth rate which jumps across a 
boundary. The half space problem was studied previously 
in a discrete version, for the (symmetrized) random per- 



mutations/PNG model 27 28] and found to also involve 



TW distributions in the limit of large system size. In order 
to exhibit full KPZ universality, it is important to solve 
the problem directly in the continuum, i.e. for the KPZ 
equation itself. Furthermore, previous approaches did not 
provide any information about the finite time behavior 
which is also universal 
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The aim of this Letter is to present a solution of the 
directed polymer problem in the continuum in presence of 
a hard wall (absorbing wall) using the Bethe ansatz (BA). 
Equivalently, we obtain the one-point height probability 
distribution for the KPZ equation on the half line x > 
with fixed large negative value of h or of — V/i (i.e. a 
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Fig. 1: Solid line: a DP with both endpoints fixed at small 
x — e with a hard wall at x — 0: the DP probability vanishes 
at the wall. Dashed line: mirror image discussed at the end. 



small contact angle) at x = 0. For simplicity we study a 
DP with both endpoints fixed - which corresponds to the 
droplet initial condition in KPZ - near the wall. We do not 
consider the case of the attractive wall although we briefly 
mention it at the end. We obtain an exact expression for 
the generating function of the moments of the DP parti- 
tion sum as a Fredholm Pfaffian, from which we extract 
the PDF of the free energy of the DP (height of KPZ) at all 
times. We then show that this PDF converges to F4, the 
Tracy Widom distribution of the largest eigenvalue of the 
Gaussian Symplectic Ensemble (GSE). The calculation is 
performed on the DP formulation, the consequences for 
the KPZ equation being detailed at the end. Our results 
are checked against numerics on a discrete DP model, both 
at high and zero temperature, thereby confirming univer- 
sality. Some consequences for extreme value statistics are 
discussed. Note that this is the first occurrence of the F4 
distribution and of the GSE within a continuum BA cal- 
culation. It is consistent with the results of 
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the discrete model and confirms that these belong to the 
same universality class than the continuum KPZ equation 
on the half space, solved here for all times. 

Directed polymer: analytical solution. We consider the 
partition function of a DP at temperature T in the contin- 
uum, i.e the sum over positive paths x(t) G R + starting 
at z(0) = y and ending at x{t) = x 

Z(x, y, t) = f X{t) ~ X Dx(r)e-*Io d T [^f + v{ x{ r)^)] ( 

Jx(0)=y 

(1) 

with initial condition Z(x, y, 0) = 5{x — y). The hard wall 
is implemented by requiring that Z(0,y,t) = Z(x,0,t) = 
0. The random potential V(x, t) is centered gaussian with 
correlator V(x, t)V(x', t) = c5[t—t')8{x—x'). The natural 
units for the continuum model are t* — 2T 5 /c 2 and x* = 



T 3 /c which allow to remove T and set c = 1 [^The time 
(i.e. polymer length) dependence is embedded in a single 
dimcnsionless parameter: 



A= (t/4t*) 1/3 



(2) 



as defined in our previous works [19 22 23 and in [20] . 
Replicating ([I]) and averaging over disorder one finds 



30 that the n-th integer moment of the DP partition 
sum can be expressed as a quantum mechanical expecta- 
tion for n particles described by the attractive Lieb-Liniger 
Hamiltonian 31 



E 
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2c 5{xi — Xj). (3) 

l<i<j<n 



in natural units (for the moment not rescaling by c, as 
in |19| ). The moments of the partition sum with both 
endpoints fixed at x can be written as: 



Z(x,x,t) n 



t 1 



-tE a 



HI 2 



(4) 



i.e. a sum over the un-normalized eigenf unctions ^ (of 
norm denoted ||/i||) of H n with energies E^. Here we used 
the fact that only symmetric (i.e. bosonic) eigenstates 
contribute. In presence of a hard wall at x — we must 
impose that ^^(xi, .., x n ) vanishes when any of the Xj 
vanishes. This case can also be solved, by a simple gen- 
eralization of the standard BA [32~l[33l . The Bethe states 



n are superpositions of plane waves 
tations P of the rapidities Xj (j 
additional summation over ±A 
for x\ < .. < x n : 



31 over all permu- 



l,..n), with here an 
The eigenfunctions read, 



^ 1 P6S„e 2 ,..£„=±l 



X^A Pl ,e2Ap 2 ,..,6 n A Pra sin(a;iAi) < 

i=2 



(6) 



A 



Ai,..,A„ 



n a 



n>£>k>l 



At — Ak Ai + Afc 



recalling that ^(a^, ..x n ) is symmetric in its arguments. 
Imposing a second boundary condition at x = L, e.g. 
also a hard wall, one gets the corresponding Bethe equa- 
which determine the possible sets of Xj. The 
and we do not repro- 



tions 
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large L limit was studied in 
duce the analysis here. The structure of the states is 
found very similar to the standard case, i.e. the general 
eigenstates are built by partitioning the n particles into a 
set of n s bound-states formed by rrij > 1 particles with 
n = Y]j—i rrij. Each bound state is a perfect string [34[ , 
i.e. a set of rapidities X^ a — kj + ~ (rrij + 1 — 2a) , where 

1 in final result performing x — > x/x*, t — ¥ t/t* and in the free 
energy F = —TlnZ, restores T dependence. 
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o=l, ...,mj labels the rapidities within the string. Such (i.e. small A) expansion of: 

eigenstates have momentum K u = rrijkj and en- — c , , „ 

8 n _ 2 p ^ J ~ 1 ] 3 z 3 = 42.99376 A (13) 
ergy = J2j=i( m j k j ~ J2 m j( m j ~ !))• Tlie difference 

with the standard case is that the states are now invariant \ nz = , — \ 3 / 2 — 0.27162097 A 3 (14) 

by a sign change of any of the momenta Xj — > —Xj, i.e. 2 y 2 

kj^-kj. J\nzjz c — ^\/27rA 3 / 2 + 0.349154645 A 3 

To simplify the problem, we restrict here to a DP with c ^ 

endpoints near the wall, i.e we define the partition sum (lnz) 3 = 0.58226188 A 3 (15) 

for x — e — (sec Fig. [IJ . ^ to 0(A 9 / 2 ) terms. The skewness of the PDF of In 2 

behaves at short time as: 

Z= lim Zv(x,x,t)/x (8) 



Then the factor involving the wave function in Q drasti- 



- , = , } lnZ } c_ /n ^ 0.079863175 A 3 / 4 (16) 



((lnz) 2 )3/ 
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cally simplifies as lim x _j. + l^/iO^i ■■x)\ 2 /x 2n = n! 2 A 2 ..A 2 . It is interesting to compare with the same results in Ref. 
The last needed factor in Q is the norm, usually not triv- 19 in the absence of the hard wall (full space) and we 



ial to obtain 35 . With some amount of heuristics we find the universal ratio of the variances: 



arrive at the following formula [36] (we now fully use the ^ — y 2 c,HS 

natural units, hence setting c = 1): p= ^"^ - ps = - - 0.076597089 A 3/2 + 0(A 3 ) (17) 



Qnz] 



11 1 1 2 lo-n. T7 o TT n r»« and of the skewness: 

\\H\\ = n\2 " J_J_Dfc 4 , m< [I U kitmitkjtmj L 3 

»=i i<i<j<n s li S lli S = 0.63689604 + 0(A 3/2 ) (18) 

m 2 fc 2 + (to + 1 - 2p) 2 /4 at small time. 

fc,m — 22m-2 11 fc2 _(_ ( m _ 2p) 2 / A We now study arbitrary time, i.e. any A, and to this 



and 



aim we define the generating function of the distribution 
P(f) of the scaled free energy \wZ = —XJ: 

f_ e -\syi 



4(fc 1 -fc 2 ) 2 + (TO 1 +TO 2 ) 2 9(a) = exp(-e-^Z) = 1 + ^ v ^ ' Z n (19) 

4(fci - k 2 ) 2 + (toi - to 2 ) 2 n_1 

from which P(/) is immediately extracted at A — > 00 



^ 4(k 1 +k 2 y + {m 1 +m 2 ) 



4(fei + fc 2 ) 2 + (mi - to 2 ) 2 lim g(s) = (?(/ + s) = Pro6(/ > -s) (20) 

We now have a starting formula for the integer moments and below we recall how it is extracted at finite A. The 

constraint X)"=i m » = ^ hi (111 can then be relaxed by 



" n'2™ s reorganizing the series according to the number of strings: 

n.=l s (m 1 ,..,ro„,) n 5 (s) = 1 + V — Z(n s , s) (21) 

f) f *3^ e (m?- m M-»i^ ff D k . m . k . m . 

- LJ - ' 2tt 4to o 15 Solvability for the generating function arises from the pfaf- 



j — 1 J l<i<j< n s 



fian identity: 



with b k , m = l[]L (ik 2 +f)- Here ("ii, ..,to„J„ stands TT n. _ TT TO j „rr ^i ~ ^ 



II £, *i>mi,fc J -,m J - - I] 1,11 v' . V J,i 



for all the partitioning of n such that Y"]".!-, to, = n with " m " J ' mj 2ifc»- X; + Xj 

° *-~*J — 1 J l<i<j<n s J — 1 J J 

to j > 1 and we used J2k ~~ ^ mjL J || which holds also (22) 

here in the large L limit. where X 2p _i = m p + 2ik p , X 2p =m p - 2ik p , p = 1, 
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This formula allows for predictions at small time. Defin- consequence of Schur's identity as used in Ref. 22 [23] to 
ing 37 z = Z/Z we obtain z 2 = z 2 — 1 as: which we refer for details. We recall that the pfaffian of an 

antisymmetric matrix A is defined as pf^4 = \J det A. Eq. 



z 2 ° = y|e 2A3 (4A 3 + 3) A 3 / 2 (erf(72A 3 / 2 ) + 1) + 2A 3 



( 22 ) allows to write the n s string partition sum as 38 
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-^A 3 / 2 + 8A 3 + 0(A 9 / 2 ) (12) mi ±J } ly 2. iik p 

v v 

X(f n T^~ m ^T t ~^ m V S T)f( 3 \ ("23") 

and, after a tedious calculation, the short time expansion " Xi + Xj > 2ns * 2ns \ > 
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Now, as in Ref. [22 23 we use the representation 
f ^ n 2S'(vi ~ Vi)e-"* Xi - v i x i = fefi and standard 
properties of the pfaffian allow to take the integral over 
the 2n s variables outside the pfaffian. After manipulations 
very similar to Ref. 22)[23 the integration and summation 
over kj, rrij can be performed, leading to: 



finds after integrations by part w.r.t. y: 

g(s) 2 = Det[7 + P K S P } = Det[J + P s/2 K P s/2 ] 
— I Ai(y + s + v, + Vj + k 2 ) 

-™ Jy>0 

Xe -i(v t -v J )k( 1 _ e iky) (30) 



K 8 (vi,Vj) 



Z(n x 



{2n s -l)\\ f = \J V]>Q 



n 



xpfp'fa-i^Wan. (24) 2->Ai{2ia)Ai{2n) 



where P x = 6(v—x) projects all Vj integrations on [x, +oo[. 
Upon using the Airy function identity 

dq 



2tt 



b)e 



iq(a-b) 



where (2n s — 1)!! = (2n s )!/(n s !2 ras ) is the number of pair- 
ing of 2n s objects, with the kernel: 



f{vi,v 2 ) 



OO 

E 

m—l 



dk {~i) m b k , 



xe 



2ir 2ik 

— m(v 1 +v 2 ) — 2ik(vi —v 2 ) 



(25) 



where we used that in the natural units t(= pr) = 4A 3 . 
Hence g(s) has now the form of a Fredholm Pfaffian and 
one shows [36] : 



g(s) = v/Dct[/ + /C] 

JC(vi,V2) = -26(v 1 )6(v 2 )d v J(v 1 ,v 2 ) 



(26) 



It is interesting that g(s) 2 is precisely the generating func- 
tion for the two independent half spaces (on each side of 
the hard wall) and that it is itself a Fredholm determinant 
(FD). Performing the rescaling Vj — > Xvj and kj — > kj/X 



leaves the result (26) unchanged with the scaled kernel: 
dk 



f(vi,v 2 ) 



x/ feA (e Ay ) 



2n 



Ai(y + s + vi+v 2 + 4fc 2 ) 



e — 2ik{vx — v 2 ) 

2ik 



(27) 



where we have used the now standard Airy trick 
/ Ai(y)e wy = e w I 3 to transform the cubic exponential in 
an exponential, together with the shift y — >• y + s + v\ + v 2 ■ 
The weight function f k (z) := ELi^^t" 2 )" 1 can be 
calculated explicitly and we find: 



2nk f 2 2 

sinn(47rfc) \ \J z y z 

- iF 2 (1; 1 - 2ik, 1 + 2ik; -1/z) 



(28) 



Eqs. (26), (27) and (28) is our main result at finite time 



for g(s). 

We now obtain the PDF of the free energy (i.e. of the 
KPZ height), first at large time, i.e. for A —> oo. Exami- 



nation of ( 28 ) leads to [36] : 

lim f k/x [e^] = -9(y)(l - cos(2fcy)) (29) 

A— ►oo 



Rescaling k — > k/2 and taking the derivative in (26) one 



we find Koiv^Vj) = -2 1 / 3 K{2 1 / 3 v l ,2 1 / 3 v j ) and upon 
rescaling of the Vi : 

Det[i - P S KP S ] , s = 2" 2/3 s 
1 



K(vi,vj) = K Ai (vi,Vj) - ~Ai(vi) 



(31) 



Ai(y- 



where Kai is the Airy Kernel KAi{vi,Vj) — J y>Q Ai(vi + 
y)Ai(v 2 + y). Our result (31) for the half space at large 
time can be compared with the full space result 18-21 
g F s{s) = Det[J - P s K Ai P s ] = F 2 (s), i.e. the GUE distri- 
bution. Hence, as compared to the two half spaces, the 



second term (projector) in (31) encodes for the effect of 
the DP configurations which in full space, cross x = at 



,(1) , 7 (2) 
J HS + ^HS 

two half spaces are statistically independent, one shows 



least once. Interestingly since Zps > Z HS + Z y ^' s and the 



from the definition ( 19 ) that 



gFs{s) < 9hs{s) 2 



(32) 



a bound valid at all times (not just for infinite A). 

We can now transform our result (31 1 into a more famil- 



iar form. Defining fC%° («i , v 2 ) — 0(vi)9(v 2 )K(vi+s 1 v 2 + s) 
we note that this Kernel can also be written as K,^° = 
B 2 S - \\B S 5){IB S \ where B s (x,y) := 9(x)At(x + y + s)%) 
one obtains via manipulations similar to Ref. 22,23, 40 
in the infinite time limit: 



g{s) 2 = Det[I - /C 5 °°] 



Det[I-Bi](l + -(l 



-{Det{I-B s 



B 2 



l<5» 



1 

2 (i| l-S 5 2 
Det{I + B s )) 2 



(33) 
(34) 

(35) 



using that = Det(I T B s )/Det(I ± B s ). Since 

Fi(s) = Det(I - B s ), F 2 (s) = Det(I - B 2 ) and using the 
definitions in J39] we obtain^] 

*a(s). 



9(s) 



2-( Fl ^ + F l( s) 



) = F 4 (s = 2- 2 ' 3 s) (36) 



Hence, to summarize, we find that for the continuum DP 
model in presence of the hard wall one can write: 



lnz = 2 2/3 A it 



(37) 



2 note that other conventions for F4 (e.g. wikipedia) differ by a 
factor \/2 
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where z — Zj Z and £t converges at large time in distribu- 
tion to the GSE Tracy Widom distribution F4. The same 
formula holds for the full space but with £ t converging at 
large time to the GUE distribution F 2 . 

We now obtain the PDF of the free energy at finite 
time. To this aim one follows the method used in 19 



It is written as a convolution, i.e. In Z = In Zq + Xuq 
is the sum of two independent random variables, where 



ItiZq has a unit Gumbel distribution (i.e. P{Zq) 



Then the PDF of u is obtained by analytical continua- 
tion p(u) = ^lmg(s)\ e xs^._ e xu +i0 +. Using (|26|, (27|) and 



( 28 1 and some complex analysis we find the tree energy 
distribution as the difference of two (complex) Fredholm 
Pfaffians (FP): 



A 



p(u) = — (y/Det[I + P K u P Q ]- ^Det[I + P K*P ] 
Ziir 



with the kernel: 



In Fig. [2] we show the convergence to the GSE TW 
distribution both for (i) T = and large polymer length 
t and (ii) at T > and large A. The agreement is very 
good. The variation for T > as a function of A is shown 
in more details in Fig. [3] where the (small) differences in 
the cumulative distributions (CDF) are shown on a larger 
scale. As in the previous figure the mean and variance of 
the numerical PDF's are adjusted to those of F4, hence 
this plot only shows variation of the shape of the PDF. 
The variance and mean are studied separately. In Fig. [4] 
we show the ratio of half space (HS) to full space (FS) vari- 
ances as a function of A. Since the two TW distributions 
have variances a Fl = 0.8131947928 and a Fl = 1.03544744, 
the ratio p should converge to the value 1.273308 at large 
time, which is apparent in the Fig. [4J up to finite i effects 
discussed there. Similarly the two TW distributions have 

skewness -y[ 2 = 0.2240842 and j[ 4 = 0.16550949 hence 
(38) 

the skewness ratio is predicted to increase from 0.636896 



K u { Vi ,Vj) = ^ J ^ J My + u + Vi+Vj+4:k 2 ) 



Sin(2( T^ )fc) [/lA(e^)+^ A (e^)] 



(39) 



^)=^^( / Wi )+ Wi )) (4o) 



-iF 2 I 1; 1 - 2ik, 1 + 2ik; - 

z 



f l k {z) = Aksmh{2kn)K ilk (2 x j -) 



(41) 
(42) 



Note that the same formula ( 38 1 with each FP replaced by 



its square, i.e. the FD, holds for the free energy associated 
to the union of the two independent half spaces. 

Numerical simulations We now perform numerical 
checks. Here we call t the (integer) polymer length. At 
high temperature, we follow [19 41|42| and define the par- 



. V(x,t) 



tition sum (PS) Z(t) = £ 7 . e T ^'■^t ' w ' ' of paths j t 
directed along the diagonal of a square lattice from (0, 0) to 
(t/2, i/2) with only (1, 0) or (0, 1) moves. We denote space 
x = (i — j)/2 and time t = i+j. An i.i.d. random number 
V(x, t) is defined at each site of the lattice (we use a unit 
centered Gaussian). The disorder averaged full space PS 

is Z = N t e pH / 2 where Nf s 



t ^ xt- — - v ) ^ S ^ ae num ber of 

paths of length i. The half space PS is obtained by sum- 
ming only on paths with x > 0, in effect equivalent to an 
absorbing wall (hard wall), with N^ s ~ 2*(2/£) 3 / 2 /^/n. 
We use the transfer matrix algorithm. It gives In z as an 
output with z = Zj Z. As was established in 



19 



41 



42 



in the high T limit at fixed A, where A = (i/(2P T )j r ^ 
for the lattice model, lnz can be directly compared - with 
no free parameter - with the analytical predictions of the 
continuum model with the same value of A, defined there 
by ([2]). In addition we also perform numerics at T = 
and compute the optimal path energy. 



at small time, Eq. (18), to 0.738604 at large time, a mod- 
erate variation. However the combination of finite size 
effects and finite sample make it hard to compute it with 
precision for small and for very large A. For A = 5 we find 
a value consistent with the above variation interval. Fi- 
nally the kurtosis are predicted to converge at large time 
to 7f 2 = 0.093448 and 7 f 4 = 0.0491952. 

Interestingly, the difference of the means /i of the 
two TW distributions (GSE and GUE) gives information 
about extreme value properties of the DP. Let us define 
p = Z HS /Z FS the probability that, in the full space prob- 
lem and for endpoints fixed at position x > the DP does 
not cross x = 0. p is defined for each disorder realization, 
with p ~ x 2 p/t for small x. Then at large time (i.e. large 
A) one has: 



lhl=2 2 / 3 A(/i^-^ 2 ) 



(43) 



where = -3.2624279 while ^ = -1.7710868. At 
small time (i.e. small A) one finds from the above results 



(and the ones in 19 
hence — lnp crosses over from 



f A 3 / 2 -0.0082964A 3 



)\np=-\ 

t 1 / 2 to - t 1 / 3 . Note 
that p is highly non self-averaging at low temperature: 
at T = it is either or 1, and a numerical study 36 
indicates that p = Prob(p = 1) decays algebraically with 
time. Computing the full distribution of p seems a hard, 
although interesting, task. 

KPZ equation: Let us now detail how our results trans- 
late in terms of the KPZ equation 



8th 



vV 2 h+^(Vh) 2 



■T){x,t) 



(44) 



where r](x,t)r](x' ,t') = R^x — x')5(t — t'), with Gaussian 
noise correlator R n (x) = D5(x). The Cole-Hopf mapping 
generally implies: 

V 



2v 



-h = InZ 



= DXi 



(45) 



Here however we must be more specific. The initial con- 
dition ([I]) corresponds to a wedge h(x, 0) = — w\x — y\ in 
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Fig. 2: Rescaled PDF of (minus) the free energy at large time, 
(i) solid line: analytical prediction 4-Fa(s). Histograms: (ii) 
in blue, ground state energy PDF (T 



2 iu with N = 10' 



samples (iii) 

-,10 



0) for a polymer t = 
in red, PDF of s = -2" 2/3 / 
for a polymer i = 2 iu at A = 6.3, with N = 10 6 samples. 
The numerical PDFs are rescaled to adjust the mean and the 
variance of F4. The variable s in all figures is called s in the 
text. 



^num(s)-^4(s) 
0.004 




Fig. 3: Convergence as a function of A: the difference between 
the numerical CDFs, F num (s), and the prediction for infinite A, 
F4(s), is plotted for A = 0.2 (in blue), 1 (in red), 3 (in yellow) 
with N = 2.10 s samples, i = 2 11 is hold fixed. For A = 0.2 a 
length i = 2 9 is also shown (dashed line) illustrating finite size 
effects. The statistical fluctuations due to finite sample N are 
visible on the figure. 



the limit w — ¥ 00, before y — > 0. Because of the hard wall 
we have that ^h(x, t) — ln(xy) + ^h(x, t) where h is not 
singular when both x and y approach zero, and the cor- 
respondence is really ^h(0,t) = InZ. Schematically the 
boundary conditions (BC) can be stated as h(0, t) = —00 
or V/i(0, t) = +00 (see more general ones below). Hence 
from J37|: 



2v 



h{0,t) = \nZ + 2 2 ^X & 



with, at large t, In Z 



,t. From 



37,38 



(46) 



D A 

„ is the same non universal constant (see discussion in 
42 ) in both HS and FS cases, the difference in In Z being 



1.5 
1.45 

1.4 
1.35 

1.3 



Q~HS 




0.1 



10 



Fig. 4: Ratio of variances p = f^j, for A varying from 0.2 
to 20. The crosses correspond to numerical data (N = 2 ■ 10 5 
samples, t — 2 , standard error estimation e = 3 • 10 -3 ). The 
dashed horizontal lines represent analytic predictions in both 
limits, § for A ->• and ^ = 1.2733 for A -> 00. The effect of 
finite t is clearly visible. It causes a small gap between these 
limits and the numerical data, which decreases as t increases. 
The solid line represents the Taylor expansion | |17[ ) globally 
rescaled to account for finite t. The right part of the graph 
shows the convergence of p at T — as a function of t: the 
upper point is t = 2 s , the lowest t = 2 10 . 



only sublinear in time, as ~ hit. 

Finally, we discuss the universality of our results. The 
BC we used here at x = is the hard wall Z = 0, which in 
the KPZ context corresponds to V/i = +00. Another stan- 
dard BC is the reflecting wall (RW) VZ = 0, i.e. V/i = 
(contact angle 7r/2). For the DP it can be achieved by con- 
sidering two symmetric half-spaces i.e. V(—x, t) = V(x, t) 
[44] . At T — there is no difference in the optimal path 
energy between the hard and reflecting wall, see Fig. [TJ 
At T > the two cases become different, since there is 
more entropy in the RW. However the longer the poly- 
mer, the closer it becomes, effectively, to the zero temper- 
ature limit. Hence we expect that although at finite time 
the two cases lead to different g(s), these become equal 
at large time. In fact all BC such that V/i > should 
converge to F4. This is consistent with the results of [28] 
translated into the T = lattice DP model (although the 
equivalent of the hard wall was not explicitly considered 
there). In the PNG model it corresponds to the absence 

We will 



of boundary source, or a weak enough source 27 



not discuss here the case of BC V/i < which leads to 
an unbinding transition. A similar transition was studied 
in the random permutation model |28| and in the PNG 
but not using the BA (see however [43]). 



11 



27 



model 

Work on that case is in progress. 

It is worth pointing out an application of our results 
to the conductance g of disordered 2D conductors deep in 



45 



the localized regime. Extending the results of Ref. 
we predict that L~ x ^hig should be distributed as F4 if 
the leads are small, separated by L, and placed near the 



p-6 



Directed polymer near a hard wall and KPZ equation in the half-space 



frontier of the sample (which occupy, say, a half space) . 
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